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Brownian Motion: Theory and Experiment 
A simple classroom measurement of the diffusion 

coefficient 
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Abstract 

. 

Brownian motion is the perpetual irregular motion exhibited by small particles immersed 
| in a fluid. Such random motion of the particles is produced by statistical fluctuations in 

■ the collisions they suffer with the molecules of the surrounding fluid. Brownian motion of 

c/3 . particles in a fluid (like milk particles in water) can be observed under a microscope. Here 

we describe a simple experimental set-up to observe Brownian motion and a method of 
determining the diffusion coefficient of the Brownian particles, based on a theory due to 
r~| ■ Smoluchowski. While looking through the microscope we focus attention on a fixed small 

Qh| volume, and record the number of particles that are trapped in that volume, at regular 

intervals of time. This gives us a time-series data, which is enough to determine the diffusion 
coefficient of the particles to a good degree of accuracy. 

O 

m 

^ ■ 1 Description of the experiment 

o 

As our system we took water containing a small number of milk particles. From our experience we 
found that freshly boiled and cooled milk, free from cream works well. A dilution of approximately 



1 drop of milk in 10 drops of water was used since this was found to be suitable for our observations. 

*-C Observations were made with a microscope (the one we used had 15x eye-piece and 45x objective). 
R" 1 ! The solution was placed between two pieces of transparencies (sealed off at the edges to prevent 
convection currents which, if not driven off, can drive things crazy) [See Fig.l]. One could see 

^ ■ Brownian motion of the milk particles if one looked at them carefully for long enough. However 
the motion is not as vigorous or spectacular as we are used to seeing in textbook pictures. 

If we focus attention on a single particle, it is expected that the particle will trace a random 
path [as in Fig. 2]. Supposing that in time t the particle is displaced from the origin to the 
position r. Then, for a Brownian particle the mean square displacement is given by < r 2 >= 2Dt 
and this defines the diffusion constant of the particle in the given fluid. A straightforward way 
of measuring the diffusion constant D would thus be to make a large number of observations on 
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the displacement of the particle in a given time t and find the mean square displacement. In our 
experiment however D turns out to be very small and so this method is not particularly easy. 
Hence we use a different method in which we estimate the diffusion constant in a rather indirect 
manner. We shall describe the theory behind this method in successive sections. 

We counted the number of particles contained in a well-defined square area in the field of 
view of the microscope at constant intervals of time. We did it 30 seconds apart and made 130 
observations. For this purpose a fine grid was made on a small piece of transparency sheet, by 
etching with a sharp edge, and attached to the eyepiece. The area actually enclosed by the grid 
was calibrated with a standard grid ( we took one that is used in pathological laboratories for 
counting blood corpuscles). The time-series data (time versus number of particles) thus obtained 
has information on both static properties, such as the distribution of number of particles in the 
observed area, and also dynamic properties such as the diffusion constant. We first look at the 
static properties and then go on to the problem of extracting the diffusion constant. 

2 Number distribution: The Poisson distribution 

We plotted the distribution of the frequency W(n) with which different numbers of particles n 
are observed within the marked region and compared it with the Poisson distribution that we 
expect on the basis of the theory which we shall briefly describe now. We consider a well-defined 
small element of volume v in a much larger volume V containing a large number of Brownian 
particles under conditions of diffusion equilibrium [See fig. 3]. The number of particles contained 
in v fluctuates with time. 

Now, if we count the number of particles contained in v at constant intervals of time r apart, 
then the frequency W(n) with which different numbers n of particles are counted will follow the 
Poisson distribution: 

W(n) = C — (1) 
n\ 

where v denotes the average number of particles contained in v. 

That the distribution will indeed have the form (1) can be proven, provided the following 
assumptions are held true in the given system: 

(i) the particles do not interfere with each other, and 

(ii) they are equally likely to be anywhere within the system. 

Surely, for our system of Brownian particles these two assumptions are really good ones. 

In Fig. 4 we show comparison of our experimental data with the Poisson Distribution, and 
find that they match pretty well, so as to justify the assumptions above. 

We would like to discuss something more before we come down to tackling the problem at hand, 
that of determining the diffusion coefficient from the time-series data obtained from observation. 

3 The Brownian particle as a random walker 

Each Brownian particle in a fluid is being kicked by fluid molecules randomly, i.e. from all sides 
with equal probability. A Brownian particle can, hence, satisfactorily be modeled as a random 
walker, i.e. a walker who moves in successive steps, each step being independent of the previous 
step. 



In the experiment we perform, we are essentially looking at Brownian motion in 2-dimensions 
(with the field of view of the microscope constituting the xy-plane). So now we consider a 
random walker in 2-dimensions. Let's look at the discrete case first. The random walker is 
moving on a 2-dimensional lattice, and can move up, down, right or left with equal probabilities. 
We assume steps of equal size a, the time for each step being r. Then if at some time t the 
walker is at the point r = (x,y), then at time t + r it can be at either of the four points 
(x + a,y),(x — a, y), (x, y + a), (x,y — a). The probability of it being at any of these four points 
is 1/4. Let Q(r, t) be the probability for the walker to be at the point r at time t. Then, 



Q(r, t) = ^[Q(x + a,y,t - r) + Q(x - a,y,t - r) + Q(x, y + a,t - r) + Q(x, y - a,t - r)]. 

Subtracting Q(r, t — r) from both sides, dividing by r, and defining D = a 2 /At, we get 

Q(r,i) - Q(r,t- t) Q(x + a,y,t-r) + Q(x - a,y,t - r) - 2Q(x, y,t-r) 

Z — ~7Z i 15. 



+ 



4r 

Q(x, y + a,t — t) + Q(x,y — a,t — r) — 2Q(x, y,t-r) 



Or, in the continuum limit, 



dQ(r,t) 
dt 



D 



d 2 Q(r,t) + d 2 Q(r,t) 



dx 2 



dy 2 



(2) 



which is the 2-dimensional diffusion equation, D being the diffusion coefficient. Solving this ( 
for initial conditions Q(r,t = 0) = 5(r), i.e. initially the particle is located at the origin) using 
Fourier transforms, we get 

1 r 2 



Q(r,t) 



AnDt 



(3) 



Thus, coming back to our system we conclude that the probability of occurrence of a Brownian 
particle at r 2 at time t when it was at r x at time t = is given by 



Q(ri;r 2 ) 



1 

AnDt 



exp 



Ifg ~ ri| 

ADt 



(4) 



4 Extraction of diffusion coefficient from the time-series 
data 

We begin by defining the probability P that a particle somewhere inside v will have come out of it 
during the time r. (N.B. As mentioned in section 3. we are into a 2-dimensional problem. Here 
we are no longer considering v as the 3-dimensional small volume, but rather as its projection 
onto the xy-plane. The symbol v mentioned henceforth is the 2-dimensional area we are focusing 
on and taking snapshots of). 

From Eq.(4) and according to the above definition of P 
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where the integration over ri is extended over all points in the interior of v while that over r 2 is 
extended over all points exterior to v. 
Alternatively, we can also write 



P 



1 


' exp 

ri,r 2 Gt) 




r 2 - ri 


2" 


AttDtv J 


4Dt 



dr 1 dr 2 



(6) 



where, now, the integration over both ri and r 2 are extended over all points inside v. (Physically, 
1 — P is the probability that a particle somewhere inside v will have still remained within v after 
a time r.) 

Now, if we take v to be a square area of side h (which can easily be ensured by using a square 
grid), then Eq.(6) leads to 



l-P 



Putting in the substitutions 



I r-h rh rh rh 

47rjDr/l 2 J J dx 1 dx 2 exp[-(x 2 - xxf/ADr] dy x dy 2 exp[-(y 2 - Vif/AD 

utions 

x h 
? = — 7= ', a = — 

we have 

1-P = 
With £i — ^ 2 — V: we have 



r 



(7) 



a 



x -= r r eX p[-(6-6) 2 ]« 2 

wir Jo Jo 



1-P = 



PCX r(i 

= dttJ d^ 

7T JO JO 



After an integration by parts, we finally have 



l-P 



2 

^Jtx Jo 



dt^e 



-e _ 



[1 - e~ a2 



(8) 



The a in Eq.(8) contains D, vide Eq.(7). Thus Eq.(8) relates P with D. 

The next question, which arises naturally, is how to find P from the experimental time-series 
data. We now show how this can be obtained from the mean square fluctuation in the difference 
in particle number in successive counts, which we will denote by ( A 2 ) av . We will need to use 
the Poisson distribution discussed in section 1. 



5 Determining P experimentally 

Let r be the experimental time interval between successive observations. Let x be the number of 
particles that remain in the volume v after time r, having started with n particles within v. 

Let y be the number of particles that enter the volume v in time r. Obviously y is independent 
of n, the initial number of particles within v, vide assumption (i). 
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We define the probability that % particles leave volume v in time r, having started with n 
particles as A\ n \ the probability that x particles remain in volume v after time r, starting with 
n particles within v as B^ n \ and the probability that y particles enter volume v in time r as E y . 

We have, with reference to the probability P defined at the beginning of section 4, 



l(n) _ ™ 



i\[n — i ! 



— P l (l-P) n -\ (9) 



which is a Bernoulli distribution. 

Now, our system is under equilibrium conditions and there is no net inward or outward flux of 
particles from volume v, i.e. the a priori probabilities for the entrance and emergence of particles 
from the volume v are equal. So E y will be the same as the probability that y particles emerge 
from the volume v in time r in general, i.e. starting with any value of n (the distribution of n 
satisfies the Poisson distribution (1) ). i.e., 

oo 

W{n) 4"'. (10) 

n=y 

From eqs. (1), (9) and (10), we eventually have 

_ e- p {v P)v 

Ey- . (11) 

which is again a Poisson distribution. Also, we have, similar to (9), 



B i U) = ZwZ ' „m pn ~ X ( X - P ) X - ( 12 ) 

and (: 

and y, which turn out as 



x\{n — x)\ 

From equations (12) and (11) we may calculate the mean and mean square deviations of x 



(x) av = n(l-P) 
((x-(x ) av f ) av = nP(l - P) 

(V)av = V P 

((y-(y )avf )av = VP. (13) 

Now, let us consider a situation in which we count n particles within the volume v in one snapshot 
and then find m particles in v in the next snapshot (i.e., after time r). If T(n, m) is the probability 
of such a transition, evidently, 

T(n,m)= £ B^E V . (14) 

x+y=m 

It can be proven, given eq. (14), that the mean and mean square deviations for the distribution 
of m are respectively the sums of those of x and y, i.e., 

(m)av = ( X ) av + { V )av \ 

( ( m - { m ) av f ) av = { ( x - ( x ) av f ) av + ( ( y - ( y ) av f ) av . 
Hence we have 

( m )av = n{\ - P) + v P ; 
( ( m - ( m ) av f ) av = nP(l - P) + v P. (15) 
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We define a quantity A n = m — n. Averaging A 2 over all possible values of m (keeping n fixed) 
using eq. (15), we get 

(A 2 ) av = P 2 [{v-nf-n] + {n + v)P (16) 

Now, we further average ( A 2 ) av over all values of n with the weight function W(n) given by eq. 
(1). We have, 

I A 2 ) = ( ( AM ) 

\ /av \ \ *— 1 n lav I av 

= 2v P. (17) 

We note that ( A 2 ) av is just the mean square fluctation in the difference in particle number 
in successive counts, and v is the mean number of particles. These are quantities that may 
easily be calculated from the experimental time-series data. Thus eq. (17) tells us how to find 
P experimentally. Finally, using eqs. (8) we can determine the diffusion coefficient from the 
experimental time series data. 

6 Experimental results 

For our experiment, P turned out as 0.326, and the diffusion coefficient D as 2.1 x 10 _12 m 2 .s _1 . 
This is almost ten times greater than the rough theoretical estimate using Einstein's formula 
[Ref. Box 1]. Given the simplicity and lack of sophistication of the experimental setup, this is a 
reasonable estimate. For obtaining better accuracy, a camera could be fitted to the eye-piece of 
the microscope in order to take snapshots at regular intervals instead of the manual procedure we 
have used. This would lead to more accurate data collection and also much more number of data 
points (we took only 130 data points). A vibration-free and tilt-free working table could also be 
used to reduce external disturbances. Finally the milk particles used in the experiment were not 
monodisperse, that is there was quite some variation in the sizes of particles. Using monodisperse 
particles would clearly be an improvement. 

7 Suggested Reading 
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Vol. 15, No. 1, p.44-52,1943 

Finally, the authors would like to acknowledge the enlightening discussions on the entire matter 
that they had with Abhishek Dhar, Supurna Sinha, Joseph Samuel and Sarasij Ray Chaudhuri 
at RRI. 



fig.1:The set up 
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[Figure caption: A better way to make observation is suggested 
We take 3 slides, cut out from OHP sheets, punch a hole in one 
them one above the other (glued with Vaseline), the punched he 
accommodate the solution.] 




gure caption: A single colloidal particle executes random rnot 



